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1.  Introduction 


Carrier  Sense  Multiple  Access  with  Collision  Detection  (CSMA/CD)  is  a 
channel  access  protocol  for  packet  broadcasting  local  area  networks .  This 
protocol  allows  user  stations  to  contend  for  time  on  a  shared  channel  by 
raultiaccessing  it  in  a  random  fashion.  In  an  attempt  to  reduce  the  frequency 
with  which  broadcast  message  packets  collide  on  the  shared  channel,  Ethernet 
(Metcalfe  and  Boggs  1976)  and  IEEE  802  (IEEE  1985)  use  a  technique  known  as  1- 
persistent  CSMA/CD  with  binarv  exponential  backoff  (Goodman  et  al .  1988).  An 
alternative  technique  proposed  for  reducing  the  frequency  of  col''’':ions  is  p- 
persistent  CSMA/CD  (see,  for  example,  Stallings  1987  and  Takagi  and  Kleinrock 
1985a);  this  protocol  is  considered  in  this  paper. 

The  dvnamic  behavior  and  performance  of  unbuffered  CSMA/CD  svstems  has 
been  studied  extensively  (e.g.,  Lam  1980.  Takagi  and  Kleinrock  1985a,  and 
Tobagi  1980).  These  studies  assumed  that  there  are  an  infinite  number  of 
network  users  and  each  user  can  have  at  most  one  packet  in  the  message  queue 
at  any  time  (called  an  infinite  source  .model).  This  simplifying  assumption  has 
been  relaxed  in  recent  studies  on  CSM^VOD  systems:  .■^postolopoulos  and 
Protonotar ios  (1986),  Goodman  et  al.  (1988),  Hastad,  Leighton  and  Rogoff 
(1987),  Park  and  Bartostynski  (1990a),  Park  and  Bartostynski  (1990b),  Takagi 
and  Kleinrock  (1985b),  and  Tasaka  (1986).  In  these  studies,  a  finite  number  of 
users  are  considered  to  have  m.essage  queues  with  finite  or  infinite  capacity 
(called  a  finite  source  .model  or  a  buffered  CS.‘L-\/CD  system). 

Am.ong  these  studies  on  buffered  CSMA/CD.  the  following  authors  considered 
p-persistent  CS.MA/CD.  Under  the  assumption  of  finite  capacity  queues, 
Apostolopoulos  and  Protonotar ios  (1986)  m.odeled  the  message  queueing  process 
as  a  two-dimensional  semi-.Markov  chain.  With  the  state  space  reduced  to  a  two- 
dimensional  space,  it  was  necessary  to  employ  an  iterative  approximation 
procedure  to  obtain  the  mean  packet  delay  in  steady  state.  The  transmission 
probability  p  was  arbitrarily  set  equal  to  the  inverse  of  the  number  of  busy- 
users  in  the  beginning  of  each  contention  period.  Thus  the  dependencv  of  the 
mean  packet  delay  on  the  control  param.eter  p  was  not  ir.'.’est  i  gated .  Takagi  and 
Kleinrock  (  1935b)  de'--eloped  a  stationary  joint  probability  generating  function 
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(PGF  hereinafter)  of  queue  lengths  for  a  system  with  rwo  users,  and  obtained 

an  explicit  formula  for  the  mean  packet  delay  based  on  that  PGF.  An  exact 

joint  PGF  of  the  queue  length  vector  for  a  system  with  more  than  two  users  was 

¥ 

first  developed  by  Park  and  Bartoszynski  (1990b).  Further,  they  obtained 
stability  conditions  analytically  based  on  the  PGF. 

In  implementing  p-persistent  CSMA/CD  the  major  problem  is  to  determine 
the  value  of  p  so  as  to  minimize  the  response  time  for  users  while  maintaining 
stability  of  the  system.  In  this  paper  an  approximation  procedure  is  developed 
to  compute  the  optimal  p  which  leads  to  a  stationary  operating  mode  and  at  the 
same  time  yields  the  minimum  mean  packet  delay  in  equilibrium.  An  optimization 
of  the  p-value  with  respect  to  the  response  time  has  not  been  studied  for  a  p- 
persistent  CS.'IA/CD  system  with  more  than  two  users,  though  it  has  been  carried 
out  for  slotted  ALOHA  by  Saadawi  and  Ephremides  (1981)  and  Sidi  and  Segall 
(1983)  . 

In  the  next  section,  the  p-persistent  CS.MA/CD  channel  access  protocol  is 
described  and  the  queueing  dynamics  under  this  protocol  is  modeled  as  a  multi¬ 
dimensional  semi-.Markov  chain.  A  stationary  PGF  of  pacV.et  backlog,  i.e.,  the 
sum  of  all  queue  lengths,  is  obtained  in  Section  :  based  on  the  joint  PGF  of 
the  queue  length  vector.  In  Section  4  the  Kolmogrov  forward  equation  is 
derived  applying  an  M/G/1  approximation  to  the  stationary  PGF  of  packet 
backlog,  and  an  iterative  procedure  is  developed  to  evaluate  the  mean  packet 
delay  in  equilibrium  and  optimize  p  with  respect  to  the  mean  delay.  A  numeric 
analysis  is  conducted  in  Section  5  and  the  analytic  results  are  compared  with 
simulation  results  demonstrating  the  accuracy  of  the  proposed  approximation 
scheme.  Section  6  contains  concluding  remarks. 

2.  A  S em i - Mar kov ian  Model  of  the  Message  Queueing  Process 

A  brief  description  of  the  p-persistent  CSMA/CD  channel  access  protocol, 
which  is  sufficient  for  the  reader  to  follow  the  subsequent  analysis,  follows. 

The  system  consists  of  a  single  seim-^er  (the  channel)  and  multiple, 
inf in i te - capac i tv  queues  of  customers  (message  packets).  The  number  of  queues, 
denoted  by  m,  can  be  as  large  as  200  in  local  area  networks.  The  channel  time 


is  slotted  with  the  slot  size  being  the  maximum  propagation  delay,  and  a  slot 
size  is  chosen  as  the  unit  of  time.  Discrete  time  index  t-O , 1 , . . .  is  used  to 
denote  slot  boundaries. 

If  the  channel  is  busy  no  users  attempt  transmission  (due  to  carrier 
sensing)  until  ongoing  transmissions  are  completed  and  the  channel  becomes 
empty  again.  Once  the  channel  becomes  empty,  every  busy  user  (i.e.,  a  user 
with  a  nonempty  queue)  persistently  attempts  at  every  c  to  transmit  a  message 
packet  from  the  top  of  its  queue.  In  a  transmission  attempt  a  user  samples  a 
random  number  from  the  uniform  distribution  over  [0,1].  If  the  number  is 
smaller  than  p.  it  starts  transmission  of  a  packet;  otherwise,  it  suppresses 
the  start  of  transmission.  Thus  whenever  the  channel  is  empty  all  busy  users 
contend  to  seize  the  channel  with  an  equal  chance.  If  only  one  among  all  busy 
users  samples  a  random  number  smaller  chan  p  at  time  c,  the  user  starts 
transmitting  a  packet  at  c  and  the  packet  is  successfully  transmitted.  If  two 
or  more  users  start  transmitting  packets  at  the  same  c,  packets  collide,  the 
collision  is  detected,  and  all  collided  transmissions  are  aborted  before 
transmissions  are  completed. 

It  is  assumed  that  in  each  slot  at  most  one  packet  arrives  at  a  queue 
(which  is  reminiscent  of  Poisson  process  assumption).  The  probability  chat  a 
packet  arrives  at  a  queue  during  a  slot  iS  denoted  by  A.  The  n'umbers  of 
packets  arriving  in  different  slots  at  the  same  user  are  independent,  and  Che 
num.bers  of  packet  arrivals  at  different  users  in  the  same  or  different  slots 
are  independent.  This  message  input  process  (called  a  geometric  arrival 
process)  implies  that  queueing  processes  at  different  users  are  symmetric; 
viz. ,  queue  lengths  at  individual  users  are  exchangeable  random  variables 
(Billingsley  1979) . 

Following  Takagi  and  Kleinrock  (1985a),  the  channel  state  dynamics  are 
modeled  as  depicted  in  Figure  1  (note  that  users  3  to  m-1  in  the  figure  are 
assum.ed  Co  be  idle  throughout  the  portion  of  time  shown  in  the  figure).  The 
channel  state  alternates  between  idle  periods  and  busy  periods.  During  an  idle 
period  all  users  have  emptv  queues,  and  during  a  busv  period  there  is  at  least 
one  busy  user.  Vich  the  assumed  m.essage  input  process,  durations  of  idle 
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periods  are  independently  and  geometrically  distributed;  hence  the  beginning 
of  an  idle  period  is  a  system  regenerative  point.  A  busy  period  is  divided 
into  a  number  of  sub(busy)periods ,  each  in  turn  consisting  of  a  transmission 
delay  followed  by  a  packet  transmission  period. 


Figure  1  about  here 


During  a  transmission  delay  no  packets  are  transmitted,  even  though  the 
channel  is  empty,  since  all  busy  users  continually  sample  random  numbers 
exceeding  p.  A  transmission  delay  terminates  at  time  t  when  one  or  more  busy 
users  sample  random  numbers  smaller  than  p  and  start  packet  transmissions.  The 
duration  of  a  successful  packet  transmission  is  proportional  to  the  length  of 
the  packet  crans.micted.  Packet  lengths  may  vary.  However,  for  analytic 
convenience,  we  use  the  expected  value  of  packet  lengths  and  denote  it  by  i. 
The  duration  of  an  unsuccessful  transmission,  which  is  the  tim.e  to  detect  and 
resolve  a  collision,  is  represented  by  2' .  In  a  baseband  system  like  Ethernet, 
it  takes  at  most  two  slots  to  detect  a  collision.  Thus  2'  is  usually  much 
shorter  than  2.  p-persistent  CSMA/CD  system  is  completely  described  by  the 
set  of  parameters,  {m , p , X , 2 , 2' ] ,  at  least  for  the  purpose  of  our  analysis. 

Throughout  the  paper  we  use  the  following  convention  in  defining 
notations;  A  boldface  capital  letter  indicates  an  .ti-vector  or  a  function  with 
vector  arguments,  except  that  E  and  P  stand  for  expectation  and  probaoility, 
respectively.  An  italic  capital  letter  with  an  underlined  superscript,  i, 
indicates  the  ith  element  of  an  ra-vector,  and  an  underlined  italic  capital 
letter  signifies  the  sum  of  all  elements  in  an  m-vector.  An  italic  capital 
letter  without  a  superscript  or  an  underscore  signifies  a  set  ot  a  function 
with  scalar  arguments. 

Denoting  the  set  of  network  users  by  I=(  1 , 2  ,  ,  .  ..ml  ,  we  define  Q^- 
iel ;  to  be  the  vector  of  queue  lengths  at  t,  and  2-”  packet 

:jcklog  at  t.  The  process  |Q^)  is  characterized  as  an  irreducible. 
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discrete  -  time  Markov  chain  whose  state  space  is  the  n-dimensional  vector  space 
iel] :  W^>0,  integer). 

A  semi-Markov  chain  model  of  the  queueing  process  is  constructed  as 
follows  (which  is  more  accurate  than  the  semi-Markov  chain  model  presented  in 
Takagi  and  Kleinrock  1985b):  The  beginning  and  the  end  of  each  idle  period  are 
regarded  as  embedded  Markov  epochs.  Within  a  busy  period,  the  boundaries  of 
subperiods  are  chosen  as  embedded  Markov  epochs.  Denoting  these  epochs  by 
k-0 , 1 , . . ,  with  k-0  at  t-0 ,  we  have  an  embedded  Markov  chain  (Q^)  which  is  a 
multidimensional  random  walk.  In  Figure  1,  the  embedded  epochs  are  shown  by  ♦, 
a^.  the  packet  backlog  at  those  epochs,  are  illustrated. 

3.  Stationary  Probability  Generating  Function  of  Packet  Backlog 

The  major  difficulty  in  analyzing  the  queueing  process  lies  in  that 

queues  at  different  users  interfere  with  each  other  through  the  shared 
channel.  This  problem  of  inceraccing  multiple  queues  cannot  be  solved  using 
classical  queueing  theory.  There  are  no  analytic  results  in  the  literature  for 
evaluating  the  mean  packet  delay  in  such  a  multidim.ensional  queueing  system, 
and  even  appro;<imatior.s  are  difficult  to  obtain  when  m>  2  (see  Kleinrock  and 
Yemini  1980)  .  In  this  paper  we  base  the  derivation  of  the  mean  packet  delay  in 
equilibrium  on  the  stationary  PGF  of  (2;^)  •  By  the  aggregation  of  (Q^)  into  a 
one-dimensional  process  (2;^)  we,  of  course,  lose  information  about  behaviors 
of  individual  queues.  The  approximation  needed  to  overcome  this  loss  of 
information  can,  however,  be  localized  in  the  delay  analysis  (as  will  be  shown 
in  Section  4)  so  as  to  generate  accurate  results  (verified  by  simulation  in 
Section  5) , 

To  describe  the  process  (2^1  we  first  define  the  following  random 
variables.  Notations  for  each  element  in  a  random  vector  and  the  sum  of  all 
elements  in  a  random  vector  are  defined  according  to  the  convention  mentioned 
earlier. 

X(n)  :  num.bers  of  packet  arrivals  at  respective  queues  over  an  n-slot  period. 
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:  numbers  of  packets  that  are  successfully  transmitted  from  respective 
queues  at  the  end  of  the  subperiod  which  begins  at  the  kth  epoch;  Note 
that  either  Z,-0  (0-vector)  or  Zr  -1  for  some  i  and  Zr  -0  for  all 
:  duration  (in  number  of  slots)  of  the  transmission  delay  in  the  sub¬ 
period  beginning  at  the  kth  epoch. 

:  duration  (in  number  of  slots)  of  the  packet  transmission  period  in  the 
subperiod  beginning  at  the  kth  epoch. 


U, 


The 


’  ‘^k+l'  ^k- 

■  { iGl ;  0^  >1 ) ,  i.e.,  the  set  of  busy  users  at  the  kth  epoch. 

{ iel :  Ct^Q  ,X-(R^)>1]  ,  i.e.,  the  set  of  users  who  are  idle  at  the  kth 
epoch,  but  become  busy  by  the  end  of  the  transmission  delay  in  the 
subperiod  beginning  at  the  kth  epoch, 
process  (2^)  can  now  be  recursively  defined  as 


2k+r  V^k  ^k’  2;,=0;  U^-  Z(k^+L^)-Z^  otherwise.  (1) 

In  this  equation  a(1)  represents  random  variable  T  1)  with  the  constraint 
Z(l)>  1.  Under  the  geomietric  input  process,  X{1)  has  a  binomial  distribution 
with  parameters  a  and  X.  Therefore,  if  (i.e.  ,  if  the  kth  epoch  is  the 

beginning  of  an  idle  period),  2;^+]^  follows  a  truncated  binomial  distribution 
given  by 

P(2;,^l-q|2;,-0)  -  [^■^]a‘?(1-A)"'^//{1-(1-,X)’^)  ,  1<  q<  m.  (2) 

Note  that  2;^.,.]^  given  2^“0  is  the  packet  backlog  at  the  start  of  a  busv  period. 

Next,  if  1,  the  random  drift  is  the  s'om  of  random  variables  X(R^)  , 

Z(L^)  and  Unfortunately,  the  latter  variables  are  stochastically 

dependent  on  each  other  in  a  complex  way.  First,  L,.  (and  hence  X(L,))  depends 
on  Z^.  Namely,  if  snd  otherwise  (i.e.,  if  Z^-1)  .  Random 

variable  Z^  in  turn  depends  on  |B^(  and  |Sj^(k^)|  since  the  success  of  a  packet 
transmission  depends  on  the  number  of  busy  users  at  the  end  of  the  transmis¬ 
sion  delay,  which  is  the  sum  of  IS-  I  and  IS.'(k.  )|.  The  latter,  13.'(Z,  )|, 
depends  on  X(R,).  Further,  R,  and  X(F-,/)  are  jointiv  dependent  on  |  b  ■.  |  ,  and 
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finally  is  uniquely  determined  given  Q^,  but  stochastically  dependent  on 

2^.  These  stochastic  dependencies  are  summarized  in  Lemma  1  below.  A  detailed 
analysis  of  these  stochastic  dependencies  and  the  resultant  conditional  PGF  of 
^k+l  found  in  Park  and  Bartoszynski  (1990b)  .  A  derivation  of 

the  latter  PGF  and  the  proof  of  Lemma  1  are  summarized  in  Appendix  I  at  the 
end  of  this  paper. 


Lemma  1 ; 


For  any  k  wi th  , 
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r  1  I 
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j2((k^)+EE 
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An  optimization  of  the  p-value  with  respect  to  the  mean  packet  delay  has 
been  studied  for  slotted  ALOHA,  as  mentioned  earlier  in  Section  1.  Compared 
with  slotted  ALOHA,  the  delay  analysis  and  optimization  of  the  p-value  is 
significantly  more  complex  for  p-persistent  CSMA/CD.  The  source  of  additional 
complexity  is  the  stochastic  dependencies  among  the  random  variables  as 
described  in  Lemma  1.  These  stochastic  dependencies  stem  from  the  carrier 
sensing  and  collision  decect ion  mechanisms  which  have  been  added  to  the 
slotted  AL0H,A  protocol  in  devising  p-persistent  CSMA/CD  (see  Szpankowski  1988 
and  Tsybakov  and  Mikhailov  1980  for  a  Markov  chain  model  of  slotted  ALOHA  with 
buffered  users).  In  this  paper  the  stochastic  dependencies  among  I, ,  Z, ,  R,, 
X(R,)  and  j3. 1  are  exacclv  incorporated  in  the  delay  analysis.  An 
approximation  is  introduced  to  simplify  the  dependency  of  \3.  |  on  2r  • 

The  process  12;,,)  is  a  random  walk  with  random  drifts  The  process  is 

uniformly  downward  bounded  meaning  that  P(2;£<'i)”0  for  all  k.  It  also 
possesses  an  important  property  called  the  bounded  homogeneity ,  which  is 
described  in  the  following  lemma  (see  Appendix  I  for  the  proof) : 


Lemma  2 :  For  any  k,  k',  and  Qj^,  such  that  I I  =  1 1  ■ 

p^l:.^.=u|2_;,.). 


Define  F,  (s)  to  be  the  PGF  of  2',-  in  view  of  (1)  we  write 

K 


Letting  W(s):-  I ’  '*’’®  ot»tain  from  Equation  (2) 


H(s)-  {  (ks+l-X)®-(l-A)'^)/|  l-d-i)'"’)  . 


By  Lemmata  1  and  2  we  may  write 


(3) 


d) 


2.  (  U.’ 

Fk^^(s\ak>0)-  eU  '"eE^s 


mini q . m) 


IS; 


0,  >0 


P(2,>0) 


)  ^  ^  P(d=ul  P(  |3..1=b|Q^=q)  P(2;^=q).  (5) 


5=1  u=-l 


Define  iJ ,  d  [  j  3 ,  |  =5) ;  =  Z*  ,  s'''P(r.,=L  I  !  3 ,  [  =5 )  .  i.e,,  the  conditional  PGF  of 

K  K  *1  A  K 

'^k  given  |S^|-5.  By  Lemma  2,  for  any  5^0,  G(s|5):-  G^(s|j3^!-5)  is  independent 
of  k  and  Qk-  By  setting  s—  =s  for  all  iel  in  Equation  (A4)  in  Appendix  I,  we 
get 

r~*  _  K  ^  •«.  ^ 

As 


C(s 


|5)=  ^ 


;=0 


[(i-p)(xs+i-,\)-a-x)J 


j 


\s-Kl-X)"^(5+j)p(l-p)^'^J'^  +(Xs  +  l-A)~^'  (l-(l-p)^'^^  -  (b+j)P(i-p)“'^-"‘; 


n-0 


(-1) 


j-n 


1-  (l-p)^'^"(l-A)'^'^'”(Xs+l-X)^''‘^ 


n 


(6) 


I'.'ote  that  index  J  in  Equation  (6)  represents  the  number  lS^(P^)i  . 

By  substituting  (6)  into  (5),  (4) -(5)  into  (3),  and  then  by  passing  k  to 
-  in  (!>'),  the  stationary  PGF  of  (2,,),  denoted  bv  F!s).  is  written  as 
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(7^ 


F(s)-  H(s)?(Q-0) 


«  mir\(q,a) 

+  \s^  2  ^7(s|b)P(|S|-b|2-q)P(2-q) 

q-1  b-1 


As  in  (7)  we  simply  drop  the  time  index  k  from  an  argument  to  indicate  its 
quantity  in  steady  state. 


4.  Delay  Analysis 

In  the  sequel,  we  write  to  denote  the  jth  partial  derivative  of 

F(s)  with  respect  to  s,  evaluated  at  s=c .  We  need  to  evaluate  F^^'^(l)  to 
obtain  the  mean  packet  backlog  at  an  em.bedded  epoch  in  steady  state,  which  in 
turn  is  used  to  evaluate  the  mean  packet  delay.  Analytic  evaluations  of  F(s) 
and  its  derivatives  are  difficult  to  obtain  mainly  because  the  conditional 
probability/  P(|Sl=bl2=c)  cannot  be  evaluated  explicitly  in  the  presence  of 
interactions  among  queues.  Therefore,  we  introduce  an  M/G/1  approximation  in 
formulating  this  probability. 

In  this  approxim.at ion,  we  assume  that  in  steady  state  the  message  queue 
at  a  user  behaves  like  a  discrete  .M/G/1  queue  independenc ly  of  other  queues. 
The  queue  at  a  user  in  the  CSMA/CD  system  inaeed  has  a  single  server  (the 
channel)  and  geometric  arrivals  (a  discrete  version  of  Poisson  arrivals)  of 
message  packets.  But  how  about  the  service  time?  The  service  time  for  a  packet 
begins  at  the  point  when  the  packet  is  promoted  to  the  top  of  its  queue.  The 
service  time  is  defined,  in  the  context  of  an  M/G/1  system,  to  consist  of:  (i) 
a  num.ber  of  sub(busy)periods ,  since  that  beginning  point,  during  each  of  which 
the  packet  is  not  transmitted  (aue  to  sampling  a  random  number  greater  than  p) 
or  is  involved  in  a  collision;  plus  (ii)  the  subperiod  in  which  the  packet  is 
finally  transmitted  with  success.  This  service  time  is  actually  affected  by 
other  busy  users'  contention  for  the  server.  Furthermore,  the  set  of  busy 
users  m.ay  change  during  the  service  time  of  a  particular  packet.  The 
ass'um.ption  of  independence  (or  "no  interference")  among  the  qi  cues  in  steady 
state  is  thus  introduced  so  that  service  tim;es  are  i.i.d.  in  steadv  state. 

With  this  as  su.mp  t  ion .  we  now  mav  write 
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P(|S|-b)- 


(8) 


where  B  is  defined  to  be  the  steady-state  probability  of  a  user  h  -i  .g  a  non¬ 
empty  queue  at  an  embedded  epoch. 

The  no  interference  assum.ption  for  the  steady-state  behavior  of  multiple 
queues  should  lead  to  a  closer  approximation  when  it  is  applied  to  p- 
persistent  CSMA/CD  than  when  applied  to  slotted  ALOHA.  This  is  because  the 
former  suffers  less  from  interference  (viz.,  collisions)  than  the  latter  due 
to  additional  controls,  carrier  sensing  and  collision  detection .  The 
approximation  in  (8)  is  motivated  by  the  results  in  Sidi  and  Segall  (1983) 
where  an  M/M/1  approximation  is  applied  to  the  slotted  ALOHA  demonstrating 
good  accuracy  over  a  wide  range  of  system  parameters.  It  should  be  noted  that 
the  approximation  is  limited  only  to  computing  ff,  and  the  rest  of  delay 
analysis  exactly  reflects  the  dynamic  behavior  of  interacting  multiple  queues. 

In''orporacing  (8)  into  (7)  and  using  the  equation  F(0)-  P(2”01p  obtain 

.~(s|!J)-  //(s)r(0l9)  +  r(s|5)l."(s|9,)-.-''i0|O  >  .  (9) 

m 

where  r(sl9):-  ^  (7(s|6)  [?]  9^(1-?)'^'^ .  (10) 

fa-1 


Taking  the  first  derivatives  of  both  sides  of  (9)  2'ields 


.(1) 


(s|«)" 


,(1) 


(s)F(0|S)  U^^'’(s|5){F(s|5)-F(0|5)  )  ./(  1-H(s|  5)  )  .  (11) 


Bv  evaluating  Equation  (11)  at  s-1  and  using  the  condition  F(l|5)-  H(l|5)-  1, 
the  equilibrium  differential  equation  (i.e.,  Kolmogrov  forward  equation  in 
equilibri’um:  see  Kleinrock  1975)  can  be  derived  as 


0  =  (l)F(0\i)  +  |5)  (  1-F(0l  H )  i  , 


(12) 
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which  can  be  rewritten  as 


F(Ol^)-  .  (13) 

Here  (1)  and  V^^'{l\6)  are  the  expected  drift  of  packet  backlog  during  an 

idle  period  and  that  during  a  sub(busy)period,  respectively.  From  (‘i)  we  get 
-  n)>/(  l-(l-X)™)  .  V^'^'^(l|9)  is  obtained  simply  by  replacing  G(slii)  by 
G^'^'^(l|5)  in  the  RHS  of  Equation  (10).  G^'^-^(l|fa)  is  the  expected  drift  of 
packet  backlog  over  a  subperiod  beginning  with  b  busy  users,  and  is  derived  in 
Appendix  II  using  Equation  (6). 

In  light  of  (12),  the  condition  '^  ( 1 1  7  )<  0  must  be  satisfied  for  the 

system  to  achieve  equilibrium.  In  Equation  (10)  we  see  that  V^‘^'^(l|$)  is  the 

marginal  expectation  of  G^'^^ClIb)  over  fa-1,.., ra,  which  in  turn  is  equivalent 

to  E(Lhl2’>0)-  Thus  the  equilibrium  equation  (12)  implies  that  a  necessary 

condition  for  stability  of  p-persistent  CSMA/CD  (i.e.,  ergodicicy  of  (Q^))  is 

aonroxi.matelv  given  by  £((/•  0-  Furthermore,  by  Fakes'  (1969)  lemma,  the 

latter  condition  is  also  sufficient  for  ergodici of  the  one -dimensional 

random  walk  1^  can  be  easily  proved  that  (Q.l  is  ergodic  if  and  only  if 

(2.  )  is  ergodic.  Therefore,  the  condition  E(2>|2-.^0)<  0  can  be  regarded  as  a 

sufficient  and  necessary  condition  for  stability.  However,  a  theoretical  proof 

of  this  assertion  is  difficult  and  has  not  been  established  in  the  literature. 

( 

The  reader  is  referred  to  Park  and  Bartoszynski  (1990b)  for  a  theoretical 
analysis  of  stability  conditions  for  p-persistent  CSMA/CD.  Some  of  the  m.ain 
results  in  that  paper  are  summarized  in  the  following  for  later  reference  in 
Section  5.  The  conjecture  below  has  been  numerically  established. 

LeiTuna  3  :  (Q^)  is  recurrent  if  '^  ( 1 1  b)“E(2^|  I  I  ”fa )  ^  0  for  all  b-l,..,m,  and 
only  if  mA  [  2-t-(  1 -p)^/l  1  -  ( 1 -p)”’l  ]  <  1. 

Conlecture :  Given  system  parameters  lm.A.£.£'}.  for  any  b  there  exist  P^(b) 
and  p..(b)  such  that  C  ^'^(l|fa)<  0  if  and  cnlv  if  p  (b}<  p<  p,.(fa). 
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Further,  (Qj.)  is  recurrent  if  p^(l)<  p<  p^{m)  since  p^(b)<  p^(b')  and 
p^(b)<  p^(b')  for  any  b>  b’ . 

Recall  that  6  is  the  probability  of  a  user  being  busy  in  steady  state, 
and  that  F(0|5)  is  the  conditional  probability  of  all  users  being  idle  in 
steady  state  given  the  value  of  8 .  It  follows  that,  under  the  assumption  of 
independent  queues,  we  must  have  1-5-  We  thus  compute  the  value  of 

8  iteratively  using  Equation  (13)  until  the  following  convergence  criterion  is 
met : 

-  { l-F(0|S)^/“) 1  <  « .  (14) 

The  function  ^+F(C | 5 1  is  monotone  increasing  in  8  and  changes  sign 

once  in  the  range  (0,1)  of  8.  Therefore,  a  bisecting  search  method  can  quickly 
find  the  value  of  8  satisfying  (14)  even  for  a  extremely  small  value  of  t. 

Let  §  be  the  value  of  8  satisfying  (14).  We  can  now  evaluate  F^'^'^d]?)  by 
applying  L'Hospital's  rule  to  (11)  and  then  replacing  the  RHS  of  Equation  (13) 
for  F(0 1 3 ) ; 


In  this  equation,  ra(/7i-l ) A^/ { 1— ( 1-A )‘^)  from  (4),  and  W^'^^(ll?)  is 

obtained  simply  by  replacing  G(s|5)  by  G'“  (l|i))  in  the  RHS  of  Equation  (10). 
A  derivation  of  G  '^(llb)  is  provided  in  Appendix  II. 

Recall  that  F^^'^(l|$)  is  the  mean  packet  backlog  observed  only  at  the 
embedded  Markov  epoch.  In  order  to  obtain  the  steady-state  mean  packet  delay 
observed  at  an  arbitrary  slot  boundary  t,  we  need  to  evaluate  the  mean  packet 
backlog  at  an  arbitrary  t  in  steady  state  (see  Takagi  and  Kleinrock  1985b  for 
a  similar  treatment).  If  t  is  in  an  idle  period,  the  backlog  is  zero.  The 
probability  that  an  arbitrary  t  is  found  in  a  busy  period  when  the  system  is 
in  steady  state  is  giv^n  by  il-F(0|5)).  Letting  2  denote  the  mean  packet 
backlog  at  an  arbitrary  t  in  steady  state,  we  have 
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2-  ii-F(Ol$) )S/r 


(16) 


where  T  is  the  expected  duration  of  a  sub(busy)period  in  steady  state,  and  S 
is  the  expected  backlog  accumulation  over  a  subperiod  in  steady  state. 

First,  to  evaluate  T,  we  note  that  the  duration  of  a  subperiod  is  i.i.d. 
for  all  subperiods  starting  with  the  same  number  of  busy  users,  due  to  the 
bounded  homogeneity  (see  Lemma  2).  Letting  S (b)  denote  the  expected  duration 
of  a  subperiod  starting  with  b  busy  users,  we  may  write  T-  (b)P( jB j-b) . 

Furthermore,  due  to  regeneracy  of  the  queueing  process  the  duration  of  a 
subperiod  starting  with  b  (>0)  busy  users  converges  a.s.  to  E(F^+L^ | | B^ | -b) 
provided  that  the  process  is  ergodic.  Consequently,  for  any  q>0  we  have 


5(b)-  E(S:,+lJ  is,  |=b)-  E< 

r 

Fi^E 

E 

^i|n 

R^, IS^(F^)  1 

|B^l-b 

1 

. 

reflecting  the  stochastic  dependencies  among  R^,  and  (see  Lemma  1). 

Define  r(s|b)  to  be  the  conditional  PGF  of  (.^,+L^)  given  |S^|-b>0.  This 
PGF  and  5(b)-  r‘’^\l|b)  are  derived  based  on  (17)  in  Appendix  III.  Using  (A9) 
in  the  .Appendix  and  Equation  (8)  we  get 


7-  ^  r'’"-^(l|b)  [■^j5^(l-5)'^'^. 
b-1 


(18) 


Turning  to  S  in  (16),  let  us  first  define  3{q)  as  the  expected  backlog 
accumulation  over  a  subperiod  beginning  with  q  outstanding  packets  so  that  S- 
z”_^^(q)P(2”<7)  •  To  derive  /3(q)  for  q>0,  we  define  •(y|q,b)  to  be  the  total 
time  (in  number  of  slots)  that  {2j.)  spends  in  state  y  during  the  subperiod 
which  starts  with  q  packets  waiting  at  b  (>1)  busy  users.  Due  to  bounded 
homogeneity,  ,J(y|q,b)  is  i.i.d.  for  all  subperiods  that  start  with  q  waiting 
packets  and  b  busy  users.  Further,  it  converges  as.  to  its  expected  value  if 
the  queueing  process  is  ergodic.  Therefore,  for  an;/  q>0  we  have 
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=0  minC  q  ,m) 

Piq)-  ^  ^  7  E?(ylq.b)  ?( lBl-bl2-q) .  (19) 

y-q  b-1 

The  process  (2j.)  is  non-decreasing  during  a  subperiod  and  the  packet  arrival 
process  is  time  homogeneous.  Consequently,  the  average  backlog  at  an  arbitrary 
slot  boundary  in  the  subperiod  which  starts  with  q  waiting  packets  and  b  busy 
users,  or  is  estimated  by  q+( .5)mXS (b) .  We  thus  can  write 

®  <«  min(q,n) 

S-  ^  ;9(q)P(2-q)-  ^  \  S  (b)  {q+(  .5)mXS  (b)  )  P  (|  3  | -b  ]  fi-q )  P  (2-c )  .  (20) 

q-1  q-1  b-1 

From  (9)  we  see  that  Z^_,qP(2-q)-  W(s  |  ^  )  ( r(  s  1  5  ) -F(0 1  5  )  )  j /ds  j  F^-^^(ll6i)  + 

q  1  1 S— I 

"^  ( 1 1  5  )  { 1 -F(  0  I  5  )  )  .  Incorporat  ing  the  latter  and  (8)  into  (20)  we  get 

n 

S-  r|^F'^-^'’(l|?)+W<'"'’(ll$){l-F(0|5)l  +  (.5).':7A  ^  F \  1 1  b)  ^  ?^  ( 1  -  5  ^ .  (21) 

b-1 

Finally,  using  Little  (1961) 's  theorem  the  mean  packet  delay  in  steady 
state,  D.  is  given  by 

D-  2/niA-  {1-F(0|$)  IS/mAT.  (22) 

The  mean  packet  delay,  D,  includes  (i)  the  elapsed  time  from  the  arrival  of  a 
packet  until  the  packet  is  promoted  to  the  top  of  the  queue,  (ii)  a  number  of 
sub(busy)periods  during  each  of  which  the  packet  is  not  transmitted  or  is 
involved  in  a  collision,  and  (iii)  the  subperiod  in  which  the  packet  is 
successfully  transmitted. 

5.  Numeric  and  Simulation  Experiments 

The  algorithm  to  compute  D  given  parameters  ( p , m, A , 2 , i ’ )  was  encoded  in 
Fortran  77;  it  first  finds  $  and  F(0l9)  satisfying  the  condition  ilA)  with 
‘’l.F-10,  and  then  compute  D  using  Equations  (15)  through  (22). 
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Computational  results  are  shown  in  Tables  1  and  2,  respectively,  for  two 
different  cases.  The  two  cases  have  the  same  values  for  m  (-50)  and  £'  (-3), 
but  different  values  for  A  and  2.  The  latter  parameters  were  set  to  1.6E-4  and 
75,  respectively,  in  Case  1  to  represent  a  high  utilization  of  the  channel 
capacity;  and  to  1.2E-4  and  25,  respectively,  in  Case  2  to  represent  a  low 
utilization.  Note  that  by  Lemma  3  an  upper  bound  on  the  channel  capacity 
(i.e.,  the  maximum  achievable  throughput  rate  in  equilibrium)  can  be  given  by 
1/2.  The  ratio  of  the  total  input  rate,  mX ,  to  1/2  is  62%  in  Case  1  and  15%  in 
Case  2 . 


Tables  1  and  2  about  here 


In  view  of  the  conjecture  in  Section  4,  the  range  [ 1 )  ,p^. (ra)  ]  of  p  that 
guarantees  stability  is  provided  in  each  table.  Among  the  steady-state 
measures  reported  in  the  tables,  the  new  notation,  S'*' ,  signifies  the  set  [b: 

0'..  The  p-'.’alues,  at  which  V*'^'^(ll§)<  '  and  S'*"  is  nonempty,  yield  a 
finite  mean  packet  delay,  though  not  satisfying  the  sufficient  condition  for 
stability  given  in  Lemjta  3.  It  should  be  noted,  however,  that  the  delay  at 
p-p,.(.73)  is  close  CO  the  minimum  delay:  in  Case  1  the  minimum  delay  was  128.6 
at  the  optimal  p-value  .16,  while  the  delay  at  p-pj^.  ( 50)- ,  08  was  140.4, 
exhibiting  a  9%  gap  from  the  minimum:  in  Case  2  the  gap  was  onlv  2%.  It  is 
also  notable  that  the  expected  n'umber  of  busy  users  in  steady  state,  or  E(|S|) 
(-  m? ) ,  was  1.5  and  1.1,  respectively,  in  Cases  1  and  2  when  p  was  set  at  the 
optimum.  This  indicates  that  there  is  little  interference  among  queues  in 
steady  state  if  p  is  optimized  with  respect  to  the  mean  packet  delay,  which 
justifies  the  M/G/1  approximation.  Indeed,  the  ratio  of  collisions  to  all 
packet  transmissions  observed  in  simulation  runs  was  about  5%  and  1%  in  Cases 
1  and  2,  respectively,  with  the  optimal  p-values. 

The  simulation  model  was  written  in  SLAM  II  (Pritsker  1984)  and  run  on  an 
IB.'i  3133.-.P  mainframe.  In  the  sim.ulation  m.odel,  we  generated  packet  arrivals 
for  each  user  so  that  the  interarrival  times  follow  an  exponential 
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distribution  with  intensity  A.  It  is  well  known  that  the  geometric  distribu¬ 
tion  of  interarrival  times  can  be  closely  approximated  by  the  exponential 
distribution  when  A  is  small  (Feller  1968).  This  exponential  approximation 
tremendously  reduced  the  computing  time.  Approximately  5,000  packet  transmis¬ 
sions  were  collected  to  estimate  the  mean  packet  delay  for  each  p-value.  The 
first  500  observations  were  discarded  to  eliminate  the  initial  transient 
period  (see  Wilson  and  Pritsker  1978,  Schruben  and  Goldsman  1985). 

The  simulation  results  are  shown  in  Figure  2  along  with  the  numerical 
approximation  results.  The  solid  line  and  dotted  line  show  approximation 
results  for  Cases  1  and  2,  respectively.  Simulation  results  for  Cases  1  and  2 
are  plotted  using  svmbols  and  ,  respectively,  and  show  good  agreement 

with  approximation  results.  The  delay  (versus  p)  curve  is  U-shaped  and  the 
minimum  delay  is  attained  in  the  neighborhood  of  the  m.aximum  p-value  that 
leads  to  equilibrium  (viz.,  near  the  bottom  right  corner  of  the  U  curve).  A 
similar,  yet  different,  pattern  of  the  delay  curv'e  was  found  for  slotted  ALOHA 
systems  by  Sidi  and  Segall  (  1983).  The  CS.MA/CD  achieves  close -minimal  mean 
packet  delays  over  a  much  wider  range  of  the  p-value  chan  the  slotted  ALOMA, 
presumably  due  to  the  additional  medium  access  ccr.Tol  mechanisms  —  carrier 
sensing  and  collision  detection. 


Figure  2  about  here 


6  .  Cone lus ion 

A  multidimensional  queueing  process  with  interactions  among  individual 
queues  arises  in  many  computer  and  communication  svstems  such  as  coupled 
processors  CSzpankowski  1988),  ALOHA  satellite  communication  (Sidi  and  Segall 
1983)  and  CSMA/CD  local  area  networks.  However,  there  are  no  analytic  results 
for  com.puting  the  mean  queueing  delay  in  such  a  svstem  with  more  than  two 
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queue s  . 


In  this  paper  we  have  shown  that  the  mean  packet  delay  in  CSMA/CD 
networks  can  be  obtained  quite  accurately  using  an  M/G/1  approximation.  The 
success  of  this  approximation  method  is  ascribed  to  the  following:  (i)  The 
message  queueing  process  has  the  bounded  homogeneity  property  (Lemjna  2).  With 
this  property,  after  introducing  the  approximate  distribution  of  the  number  of 
busy  users  (Equation  (8)),  we  could  simplify  the  stationary  PGF  of  the  packet 
backlog  (Equation  (7))  into  a  product  form  (Equation  (9)).  (Note  that  in 
Equation  (9),  S*_^P(G-q)  in  Equation  (7)  is  factored  out  as  ( T( s | 9 ) -F (0 | 9 ) } . ) 
(ii)  The  delay  analysis  incorporates  an  exact  PGF  for  the  random  drift  of  the 
packet  backlog  over  a  period  between  two  consecutive  embedded  Markov  epochs 
(Equations  (2),  (6)  and  Lemma  1).  (iii)  Finally,  the  interference  among  queues 
is  actually  minimal  in  steady  state  if  the  control  parameter  p  is  optimized. 

A  practical  and  easy  way  to  determine  the  p-value  would  be  to  set  it 
equal  to  p^^.(m);  this  guarantees  stability  (by  Lemma  3  and  Conjecture)  and  at 
the  sam.e  time  yields  a  close-minimal  mean  packet  delay  in  equilibri'um  (as 
numerically  dem.onstrated  in  Section  5).  According  to  Conjecture,  for  given 
system  parameters  p,.(m)  is  obtained  !".■  solving  f  1  |.m)-  0  for  p. 

Substituting  a  for  5  in  Equation  (,A5)  in  App-:'.dix  II,  we  obtain 


G'*-’(l|.m)- 


1  -  ( 1  -  o  ) ' 


i(l-p)•■^-^(mii-l)p(l-p)'^■‘-^ii’  (1-  (l-p)‘"-.mp(l-p)'’^'^)l 


(23) 


Rewriting  (23)  we  see  that  G^'^'^(l|n)-  0  is  equivalent  t 


/(q):-  -Ciq‘’^-t-Ciq-’^' ^+C3-  0, 


(24) 


where  q-l-p,  Ci-  c^  +  A ( i '  - 1 )  ,  Cj-  mA(i-£')-l  and  €3=  XI'.  The  equation  /(q)  has 
the  global  minim.um  at  q  -(m-l)c2/n7c,  .  It  has  exactly  two  solutions  (qi.q^) 

"if 

such  that  0<  q^K  q  <  q2<  1  provided  that  mXl.<  1  (which  is  a  necessary 
condition  for  stability  by  Lemjna  3)  and  f(q  )<  0.  A  bisecting  search  over  the 
range  (O.q'^)  of  q  can  quicklv  find  q,.  and  then  we  have  p.^.Lm)-  1-q. 
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Our  approach  --  viz.,  incorporating  an  exact  transient  analysis  of 
multiple  interacting  queues  into  an  M/G/1  approximation  for  evaluating  average 
delay  in  steady  state  --  can  be  employed  to  analyze  other  CSMA/CD  protocols. 

It  would  be  Interesting  to  attempt  applying  this  approach  to  the  analysis  of 
new  variants  of  CSMA/CD  that  are  designed  for  high-speed  broadcast  bus 
networks  (e.g.,  p ^  persistent  protocol  by  Mukherjee  and  Meditch  (1988),  LCSHA- 
CD  by  Maxemchuk  (1988),  and  a  modified  CSMA/CD  by  Lin  and  Sousa  (1990)). 
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Appendix  I.  The  Conditional  Joint  ?GF  of  Given 

The  joint  probability  of  arid  X(r)-[ ,  iel :  conditional  on 

is  given  by 


?(R^-r,X(r)-[x.,ieI]\q^)-  -[l-(l-p) 


(1-1)4  {  M  ) 

'(r)  J  ^ie3kir)^i^ 


f  ,  {r-w  '\  x.-l  r-w-x  -t-l-v 

Y  (1-A)'^‘4(1-p)^'’*'|^^  ^  (1-1)  "  j. 


In  this  equation,  w  represents  the  slot  in  which  the  first  packet  since  the 
beginning  of  the  kth  subperiod  arrives  at  user  i  (who  had  no  waiting  packets 
in  the  beginning  of  the  subperiod).  It  should  be  noted  that  B.^(r)  is  uniquely 
determined  given  and  X(r).  and  is  uniquely  determined  given  Q^.  Equation 
(Al)  thus  shows  that  the  joint  distribution  of  and  X(i?^)  depends  only  on 
3^,  apart  from  the  system  parameters  p.  1.  and  I.  That  is,  the  joint 
distribution  is  identical  for  anv  k  and  k'  with  3  =3,,,  even  if  /. 

The  probability  of  having  a  collision  in  the  subperiod  beginning  at  the 
kth  epoch,  given  is 


|5.  |  +  |B’ (r) I 


|3^|HSr(r)|-l. 


r  |5;,i  +  is;(r)h 

l(l-(l-p) 


The  probability  of  having  a  successful  transmission  in  the  subperiod  beginning 
at  the  kth  epoch,  conditional  on  is  simply  given  by  P(Z^i^0  ] , B^( r) )-  1- 

P  (Z. -0  I  3 .  ,  B ,' ( r ) )  .  As  such,  given  P.-r,  Z,  depends  only  on  B,  and  31(r).  Recall 
that  B'(R,)  is  uniquely  determined  by  3,  and  X(R,)  and  that  R,  and  XiR,)  are 

K  K  K  K  K  K 

jointl'.'  dependent  on  S^.  Therefore,  Z^,  is  transitively  dependent  only  on  3^. 
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on 


Since  X(L^)  depends  only  on  Z^.  it  follows  that  the  random  drift  depends 
only  through  We  thus  have  proved  both  Lemmata  1  and  2. 

Let  f ( s  1  ,  where  s-[s,.  ,iel],  be  the  conditional  PGF  of  given 

Q^»<0 ,  and  G(s|S);-  G^(sl3^-3)  be  the  conditional  PGF  of  given  which  is 

independent  of  k  and  (by  Lemma  2) .  Then  we  have 


k+l 


(s|q^^0)= 


of.  {  I'f. 

1 

E 

1  1  1 1 
[i€l  [i&I 

’4 

Q.  ^0 

P(Q 


v''®  )  1  US'-  r 


G  ( s  I  3 ) 


b=l  VSCI;  t31=h 


F^(s:s.-0  Vi^)-  ^  F^(s:s.  ,-0,s.-0  Vi^)]  .  (A3) 

i  '£3 


Using  (Al)  and  (A2),  function  G(sj3)  is  evaluated  as: 


G(sl3)=  e] 


r 

s  .  E! 


15.  E  5 


F..X(Fp 


1 


X  ; 
5  .  * 


-=»o  X .  =0  Vi€r 

i 


Z..-0  ViGi 


iel 


Z,^' z  .  . 
<  1 


iel 


■z :] 

S:  HX 


X  P 


Z. =[z  .  , i£l’ 

K  "  1 


'' 

F.  =r,X(r)-;x,.  ,  iei;  .3.  =3  P  F.  =r ,  X  ( r )  -  ;  x  ,  isl  ] 


H- 


n-  b 


■I  }  ( 

j-0  VS'Cl-3:  |B'  |=j 


I  I 

'  =3  ' 


t1 


(1-p)  (is^  +  l-  X)  -  i^l-X) 


i^3uB 


+  $(s)(l-(l-p)^^J-(5+j)p(l-p)^^J‘^)|  X 


X  2  ^  (-1) 

.■;-0  VAC3'.lAl=n 


j-n 


l-(l-p)“‘^'’(l-X)-'^'^''^  M  (As.+l-X) 

i&Su'A 


-1 


(A4) 


h  O 


where  n(s)-  ||  (As.  +  l-A)"^  and  ^Cs)-  ||  (as.+I-A)"^  are  the  PGF's  of  the 
iel  ^  iel  ^ 

duration  of  a  successful  and  an  unsuccessful  transmission  period, 
respectively.  A  detailed  derivation  of  Equation  (A4)  can  be  found  in  Park  and 
Bartoszynski  (1990b) . 


Appendix  II.  Derivatives  of  the  PGF  CCs|b’)  Evaluated  at  s-1 
Taking  the  first  and  second  derivatives  of  G(s|b),  which  is  given  in 
Equation  (6),  with  respect  to  s  and  evaluating  the  derivatives  at  s-1  we 
obtain 


V  .m-b)  (  A  r 

( 1 1  b )  -  \  -  0 

.4.  ;  i  -  n 


J  J  u- 


(l-d-p)-^  -')jp(l-A)/(p-A)  +  (Anji-l)a  + 


+  Xmr  !  1-  (1-p)' 


G^-^dlb) 


- k.'a(  A‘.-::luJid)-2A.mi+2)  +  i 

A-pJ  i- 


-.77i'  (mi'  -1)  d-  (l-p)^'^-^-Qll  + 


+  2  (A.m2-1)Q  +  Ami'  (1-  (1 


-p)^^-^'  -a) 


A.^  +  Jp(l-A)v/(p-A)  -i- 


f  !  1  -  (1 -p)  )  A^(2ir+0')  +  2JpA  (1  -  A  )  d/ (  p  -  A  "I  f 


+  jd-d-p)^'"-''lv  j-1  +  2jA(l-p)/(p-A)  +  (J  +  l)A2(l-p)2/(p.;^)2l  _ 


where  a:-  (b+j  )p(l ‘  ^  1- ( 1 -p)^'*'"d  -  A)'^' ^ 


■  7  n 


J  ( -l)J'”(b+.n)  (b+.nd)  d-^g),^-; 
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II  r\^8 


-  ill!) 


Appendix  III.  Derivation  of  the  PGF  r(slb’)  and  Its  First  Derivative 
In  this  Appendix,  a,  0  and  ^  defined  in  Appendix  II  are  used  again.  Based 
on  Equation  (17)  we  can  write 

r(sli))-  Ejs'^^  E  E  |S^(i?^)  ||  IS^I-bj- 

m-hf  bl  ^  r  L  1 

s'"  ^  I  ^  E  s  p[z^=2||Bi(r)l-J.|B^l-b]p[i?^-r,|B'(r)|-j||B^l-b]- 

^  0  j-0''  J  Jz-0  L 


®  m-b  t 


r-0  j-0^  ^ 


Evaluating  the  suirjnation  over  w  and  then  the  infinite  series  over  r,  we  obtain 


r(s|b)- 


^  —  pQ  +  s'^'{l-(l-p)^'^-^'-a)]{  )  P](-l)^’‘''/{l-s(l-^))},(.A8) 

I  J  U-pJ  ^  ^  nin^‘ 


r^-^^dlb)-  ^  I"’/'  — j  |v-[icr  +  £'  -a)  +  (l-Cl-p)^"^^)  X 


X  j  )  d  (-l)^'"(l-^)/;3' 
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Figure  1 .  Channel  States  and  Queueing  Process  in  a  Slotted  p-Persistent  CSMA/CD  System 
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Legend:  A  generation  of  a  new  packet 

X  transmission  attempt  sampling  a  random  number  greater  than  or  equal  to  p 
p  transmission  attempt  sampling  a  random  number  smaller  than  p 
n  total  packet  backlog  in  the  system  at  the  kth  embedded  Makov  epoch 
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